Melting behavior of large disordered sodium clusters. 
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The melting-like transition in disordered sodium clusters Najv, with N=92 and 142 is studied 
by using a first-principles constant-energy molecular dynamics simulation method. Nai42, whose 
atoms are distributed in two (surface and inner) main shells with different radial distances to the 
center of mass of the cluster, melts in two steps: the first one, at « 130 K, is characterized by a 
high intrashell mobility of the atoms, and the second, homogeneous melting, at ~ 270 K, involves 
diffusive motion of all the atoms across the whole cluster volume (both intrashell and intershell 
displacements are allowed). On the contrary, the melting of Na92 proceeds gradually over a very 
wide temperature interval, without any abrupt step visible in the thermal or structural melting 
indicators. The occurrence of well defined steps in the melting transition is then shown to be related 
to the existence of a distribution of the atoms in shells. Thereby we propose a necessary condition 
for a cluster to be considered rigorously amorphouslike (totally disordered), namely that there are 
no space regions of the cluster where the local value of the atomic density is considerably reduced. 
Nag2 is the only cluster from the two considered that verifies this condition, so its thermal behavior 
can be considered as representative of that expected for amorphous clusters. Nai42, on the other 
hand, has a discernible atomic shell structure and should be considered instead as just partially 
disordered. The thermal behavior of these two clusters is also compared to that of icosahedral 
(totally ordered) sodium clusters of the same sizes. 

PACS numbers: 36.40.Ei 64.70.Dv 



I. INTRODUCTION 

Cluster mftLtpg is a topic of current theoreticalEFB and 
experimentaflili interest, motivated by the observation 
of several size-dependent properties which have no ana- 
log in the bulk phase. Between those properties, we can 
mention the following: (i) the melting like transition does 
not occur at a well defined temperature as in the solid, 
but spreads over a finite temperature interval that widens 
as the cluster size decreases. The lower end of that in- 
terval defines the freezing temperature below which 
the cluster is completely solidlike, their constituent atoms 
just vibrating about their equilibrium positions. The up- 
per part of the interval defines the melting temperature 
above which all the atoms can diffuse across the clus- 
ter and the "liquid" phase is completely established. Be- 
tween those two-temperatures the cluster is not fully solid 
nor fully liquidJl3 It is in that transition region where the 
cluster-specific behavior emerges: (ia) Premelting effects 
like partial meltinj|_pf the cluster (the most usual case 
is surface melting)!! 3 ] or structural isomesizations upon 
heating, which lead to a melting in steps.til (ib) The dy- 
namic coexistence regime, where the cluster can flucttu 
ate in time between being completely solid or liquid .LSI 
(ii) Strong nonmonotonic variations of the melting tem- 
perature with size have been found in recent experiments 
on sodium clusters. El The maxima in the melting temper- 
ature are not in exact correspondence with cither elec- 
tronic or atomic shell closing numbers, but bracketed by 
the two, suggesting that both effects are relevant to the 
melting process. It is important to note that the values 
of Tf and T m as defined above are not yet amenable 



to the experiments, and that the experimental melting 
temperature js-.somewhere between these two values. 

PreviouslyBQ we have reported density functional 
orbital-frcc molecular dynamics (OFMD) simulations 
of the melting process in sodium clusters N«wv, with 
N=8,20, 55,92, and 142. The OFMD technique^ is com- 
pletely analogous to the method devised by Car and 
Parrinello (CPl_to perform dynamic simulations at an 
ab initio level,EZl but -the electron density is taken as 
the dynamic pSiariableJiS as opposed to the Kohn-Sham 
(KS) orbitalsEJ in the original CP method. T-his-tech- 
nique has, been already used both in solid stateEjO and 
clusterBE3'E3 physics. In contrast to simulations which 
use empirical interatomic potentials, the detailed elec- 
tronic structure and the electronic contribution to the 
energy and the forces on the ions are recalculated effi- 
ciently every atomic time-step. The main advantage over 
KS-based methods is that the computational effort to 
update the electronic system increases linearly with the 
cluster size N, in contrast to the N 3 scaling of orbital- 
based methods. Indeed, these were the first molecular 
dynamics simulations of large clusters as Nag2 and Nai42 
that included an explicit treatment of the electronic de- 
grees of freedom. 

A very important issue in the simulations of cluster 
melting is the election of the low-temperature isomer to 
be heated. A good knowledge of the ground state struc- 
ture (global minimum) is required, as the details of thci 
melting transition are known to be isomer-depcndcnt.Efl 
But the problem of performing a realistic global optimiza- 
tion search is exponentially difficult as size increases, so 
finding the global minima of Nag2 and Nai42 becomes im- 
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practical. In our previous work@ we directly started-irom 
icosahedraJ-isomers, as there is some experimentala and 
theoreticalca indications that suggest icosahedral packing 
in sodium clusters, and found a good agreement with the 
experimental results of Haberland's group.0 However, we 
were unable to find those icosahedral structures by an un- 
constrained search method as simulated annealing, which 
always led to disordered isomers both for Na92 and Nai42 . 
Although the icosahedral structures are more stable in all 
the cases, the energy difference between both isomers is 
approximately 0.02 eV/atom, which is very small. Amor- 
phouslike structures have been found recently to be the 
ground rotate isomers of gold clusters for a number of 
sizeSjEa-EH and pair potential calculations performed by 
Doye and Wales predict that thp-amorphous state is fa- 
vored by long potential rangesa The specific features 
of those structures are little or no spatial symmetry and 
a pair distribution function typical of glasses. Besides 
that, one usually finds a large number of amorphouslikc 
isomers nearly degenerate in energy, which suggests that 
they occupy a substantial fraction of the phase space 
available to the cluster. Both the energy proximity to 
the more stable icosahedral isomers and the large en- 
tropy associated with the amorphous part of the phase 
space make plausible that amorphous isomers could be 
present in the experimental cluster beams, so their ther- 
mal properties deserve specific investigation. Apart from 
this, the study of melting in amorphouslikc clusters is 
intrinsically interesting from a theoretical point of view. 
Thus, the goals of this work are to study the mecha- 
nisms by which the melting-like transition proceeds in 
two disordered isomers of Nag2 and Nai42 , to study the 
influence on the melting behavior of the disorder degree, 
and to make a meaningful comparison with the melting 
behavior of the corresponding icosahedral isomers. In the 
next section we briefly present some technical details of 
the method. The results are presented and discussed in 
section III and, finally, section IV summarizes our main 
conclusions. 



II. THEORY 

The orbital-free molecular dvaamics method is a Car- 
Parrinello total energy schemeO which uses an explicit 
kinetic-energy functional of the electron density, and has 
the electron density as the dynamical variable, as op- 
posed to the KS single particle wavefunctions. The 
main features of the energy functional and the calcula- 
tion schejae, have been described at length in previous 
work,li§E3Em ancLdetails of our method are as described 
by Aguado et al.B In brief, the electronic kinetic energy 
functional of the electron density, n(r), corresponds to 
the gradient expansion, arxamd the homogeneous limit 
through second ordertSc3~Ell 



where the first term is the Thomas-Fermi functional 
(Hartree atomic units have been used) 
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and the second is the lowest order gradient correction, 
where T w , the von Weizsacker term, is given by 
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The local d en s it y approximation is used for exchange and 
correlation.E2lE2l In the external field acting on the elec- 
trons, V ext (r) = ^2 n v{r — Rn), we take v to be the lo- 
cal pseudopotential of Fiolhais et aZ.,Eil which reproduces 
well the properties of bulk sodium and has been.-shown 
to have good transferability to sodium clusters. L3 The 
cluster is placed in a unit cell of a cubic superlattice, and 
the set of plane waves periodic in the superlattice is used 
as a basis set to expand the valence density. Following 
Car and ParrinelloJlJ the coefficients of that expansion 
are regarded as generalized coordinates of a set of ficti- 
tious classical particles, and the corresponding Lagrange 
equations of motion for the ions and the electron density 
distribution are solved as described in Ref. ^|. 

The calculations used a supercell of edge 71 a.u. and 
the energy cut-off in the plane wave expansion of the den- 
sity was 8 Ryd_In all cases, a 64x64x64 grid was used. 
Previous tests&Q indicate that the cut-offs used give good 
convergence of bond lengths and binding energies. The 
fictitious mass associated with the electron density coef- 
ficients ranged between 1.0x10 s and 3.3x10 s a.u., and 
the equations of motion were integrated using the Verlet 
algorithmic for both electrons and ions with a time step 
ranging from At = 0.73 x 10~ 15 sec. for the simulations 
performed at the lowest temperatures, to At = 0.34 x 
10~ 15 sec. for those at the highest ones. These choices 
resulted in a conservation of the total energy better than 
0.1 %. 

Several molecular dynamics simulation runs at differ- 
ent constant energies were performed in order to obtain 
the caloric curve for each cluster. The initial positions 
of the atoms for the first run were taken by slightly de- 
forming the equilibrium low temperature geometry of the 
cluster. The final configuration of each run served as the 
starting geometry for the next run at a different energy. 
The initial velocities for every new run were obtained by 
scaling the final velocities of the preceding run. The total 
simulation times varied between 8 and 18 ps for each run 
at constant energy. 

A number of indicators to locate the melting-like tran- 
sition were employed. Namely, the specific heat defined 
byB 
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where N is the number of atoms and <>t indicates the 
average along a trajectory; the time evolution of the dis- 
tance between each atom and the center of mass of the 
cluster 

n(t) =| Ri(t)-R cm {t) |; (5) 

the average over a whole dynamical trajectory of the ra- 
dial atomic distribution function g(r), defined by 

dN at = g{r)dr (6) 

where dN at (r) is the number of atoms at distances from 
the center of mass between r and r + dr; and finally the 
diffusion coefficient 



which is obtained from the long time behavior 
of the mean square displacement < r 2 (t) >= 

WK; YTjli Ei=i[-Ri(*Oj + t) - Ri(toj)} 2 , where n t is the 
number of time origins, t j , considered along a trajectory. 

III. RESULTS 

The most stable disordered structures of Nag2 and 
Nai42 that we have found with the simulated annealing 
technique are shown in Fig. 1, together with the relaxed 
icosahedral isomers obtained as explained in our previous 
work.Q The dynamical annealing runalJ were started from 
high-temperature liquid isomers thermalized during 30 ps 
at 600 K. The cooling strategy was to reduce the inter- 
nal cluster temperature by a factor of 0.9999 of its actual 
value each twelve atomic time steps. With the chosen 
time step of 0.34 x 10~ 15 sec. the temperature reduction 
is applied each four femtoseconds. A first important dif- 
ference between the icosahedral and disordered isomers 
is that the former ones have a smoother surface. Besides 
that, no apparent spatial symmetry is observed in the 
disordered isomers. In Fig. 2 we show the short-time av- 
erages (sta) of the distances between each atom and the 
center of mass of the cluster for both isomers of Nai42, 
obtained from an OFMD run at a low temperature. The 
values of < 7j(t) > s ta are almost independent of time. 
The clusters are solid, the atoms merely vibrating around 
their equilibrium positions. Curve crossings are due to 
oscillatory motion and slight structural relaxations rather 
than diffusive motion. The difference between Figs. 2a 
and 2b is due to structure. For the icosahedral isomer 
in Fig. 2a one can distinguish quasidegenerate groups 
which are characteristic of the symmetry: one line near 
the center of mass of the cluster identifies the central 
atom, whose position does not exactly coincide with the 
center of mass because of the location of the five sur- 
face vacancies (147 atoms are needed to form a complete 
three-shell icosahedron) ; 12 lines correspond to the first 
icosahedral shell; another 42 complete the second shell, 



within which we can distinguish the 12 vertex atoms from 
the rest because their distances to the center of mass are 
slightly larger; finally, 82 lines describe the external shell, 
where again we can distinguish the 7 vertex atoms from 
the rest. In contrast, the lines for the disordered isomer 
in Fig. 2b are quite dispersed. Nevertheless, there is a 
narrow interval where the ionic density is very low, that 
serves to define separate surface and inner atomic shells. 
The case of Na92 (not shown) is similar, but in that case 
the structure is more uniformously amorphous, and there 
is no way to distinguish between surface and inner shells. 
We will see below thatj-this small difference is very im- 
portant. Soler at alrB'LJ have compared the structures 
of icosahedral and amorphous gold clusters, with similar 
results: while the atoms in the icosahedral isomers are 
clearly distributed in atomic shells with different radial 
distances to the center of mass of the cluster, in the case 
of amorphous clusters there are "atomic transfers" be- 
tween shells that blur the atomic shell structure. In the 
case of gold clusters, the amorphous isomers turn out 
to be-the minimum energy structures for a number of 
sizesBj In the case of Nag2 and Nai42, the icosahedral 
isomers are more stable than the lowest energy disor- 
dered isomers found (0.017 eV/atem and 0.020 eV/atom 
for Na92 and Nai42 respectively) .0 

For each cluster we have calculated the internal tem- 
perature, defined as the average of the ionic kinetic 
energy ,EZl as a function of the total energy- the so-called 
caloric curve. A thermal phase transition is indicated in 
the caloric curve by a change of slope, the slope being 
the specific heat; the height of the step gives an estimate 
of the latent heat of fusion. However, melting processes 
are more easily recognized as peaks in the specific heat as 
a function of temperature, that has been calculated di- 
rectly from eq. (4) . The fact that the specific heat peaks 
occur at the same temperatures as the slope changes of 
the caloric curve (see curves below) gives us confidence 
on the convergence of our results, as those two quantities 
have been obtained in different ways. 

The specific heat curves for Nai42 (fig. 3) display two 
main peaks, suggestive of a two-step melting process. For 
the disordered cluster the two peaks are well separated 
in temperature, Tf s m 130 K and T* 5 » 270 K, whereas 
they are much closer together for the icosahedral cluster, 
T l s c ° « 240 K and « 270 K, so close that only one 
slope change in the caloric curve can be distinguished in 
this case. The results suggest that the melting transi- 
tion starts at a temperature T s and finishes at T m with 
the difference in the melting of the two isomers being the 
smalkr T s value of the disordered isomer. In our previous 
workul we showed that for the icosahedral cluster those 
two steps are due to surface melting and homogeneous 
melting, respectively. Here we show that a similar expla- 
nation is valid for the disordered Nai42 isomer. At T=160 
K, a temperature between T* s and T* s , the structure of 
the disordered cluster is more fluid than at low temper- 
ature (fig. 2). Fig. 4 indicates that the spherical atomic 
shells have separately melted, atoms undergoing diffu- 
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sive motions mainly across a given shell, as seen in the 
bold lines, without an appreciable intershell diffusion (al- 
though some occasional intershell displacement has been 
observed) . The larger spread of the upper bold line in- 
dicates that diffusion is appreciably faster in the surface 
shell. Thus the transition at 130 K can be identified with 
intrashell melting. Why it occurs at a rather low tem- 
perature is associated with the large spread in the radial 
distances of the atoms in each shell. The atomic shells 
are structurally disordered at low temperature, although 
kinetically frozen, like in a typical glass, so inducing diffu- 
sion in the shells of that cluster is rather easy and occurs 
at moderate temperatures. The surface melting stage 
does not develop in the icosahedral isomer until a tem- 
perature of T* co w 240 K is reached. At this temperature, 
the time evolution of the surface atomic distances to the 
cluster center becomes very similar to those of the dis- 
ordered Nai42 isomer at 160 K.Q Inducing diffusion in 
the surface of the icosahedral isomer requires a higher 
temperature because of the higher structural order of its 
surface. Once the surface of both isomers has melted, ho- 
mogeneous melting occur at the same temperature, « 270 
K, in veijy good agreement with the experimental value 
of 280 KJj From that temperature on, the liquid phase is 
completely established (all atoms can diffuse across the 
whole cluster volume)Q and differences between both iso- 
mers have no sense anymore. 

The thermal behavior of both isomers is not so differ- 
ent. The radial atomic density distribution of the disor- 
dered isomer at 160 K (Fig. 5) shows a smoothed shape 
with respect to that found at low T, but a distribution 
of the atoms in separate surface and inner atomic shells 
can be clearly distinguished. In fact, the intermediate 
temperature distribution is similar to that found for.-jthc 
icosahedral isomer in the same temperature region,Q as 
the atomic shells are equally distinguished. The small 
gap in the ionic density (Figs. 2 and 5) of the disordered 
isomer at low temperature drives the cluster towards a 
well defined shell structure upon heating. We conclude 
that, despite the high orientational disorder in both sur- 
face and inner shells, the cluster can not be considered 
completely amorphous, as Fig. 5 at intermediate tem- 
perature shows some radial atomic ordering. There are 
still more similarities. The solidlike phase of the icosahe- 
dral isomer disappears as soon as a temperature of 130 
K is reached, even though no peak in the specific heat 
is detected: there are isomerization transitions between 
different permutational isomers which preserve the icosa- 
hedral structured These isomerizations involve the dis- 
placement of the five surface vacancies across the outer 
shell. Thus, both isomers leave the solidlike phase at w 
130 K, the only difference being that one has direct access 
to the intrashell melting stage while the other entexs first 
an isomerization regime. Calvo and SpiegelmannQ have 
related the appareance of specific heat peaks to sudden 
openings of the available phase space. In the isomer- 
ization regime the icosahedral cluster has access just to 
a limited number of symmetry-equivalent isomers, while 



the phase structure of the disordered cluster opens sud- 
denly to include a very large number of isomers, as all 
the possible position interchanges between two atoms of 
a given shell are allowed. Thus, a specific heat peak ap- 
pears at T«T30 K for the disordered isomer, but not for 
the icosahedral isomer. Any atomic shell distribution in 
the time average of g(r) disappears completely when ho- 
mogeneous melting occur (Fig. 5). 

The results for Nag2 are shown in fig. 6. Two-step 
melting is again observed in the icosahedral isomer, with 
a small prepeak in the specific heat at T* co as 130 K 
and a large peak, corresponding to homogeneous melt- 
ing, at T^° « 240 K. In this case T\ co and T^° are well 
separated. T* co is in the range where the isomerization 
processes in icosahedral Nai42 set in and the intrashell 
melting stage in disordered Nai42 develops. The larger 
number of vacancies in the surface shell of icosahedral 
Nag2, as compared to icosahedral Nai42, allows for true 
surface diffusion and these processes give rise to a dis- 
tinct peak in the specific heatE The disordered isomer 
melts gradually over a wide temperature interval, and 
no prominent specific heat peaks nor important slope 
changes in the caloric curve are detected. Ercolessi et 
al.cB have also found a melting process without a latent 
heat of fusion for amorphous gold clusters with less than 
~ 90 atoms. In Fig. 7 we show the radial ionic density 
distribution of disordered Nag2 at several representative 
temperatures. At a temperature as low as 70 K there is 
no discernible atomic shell structure. The g(r) function 
for a higher temperature where the surface of the icosa- 
hedral isomer has already melted is not very different. At 
a temperature where the icosahedral isomer is liquid the 
only appreciable change in g(r) is due to the thermal ex- 
pansion of the cluster. The structure of cold disordered 
Nag2 is both radially and orientationally disordered. The 
structure is kinetically frozen, but there seems to be no 
barriers impeding the exploration of the liquid part of the 
phase space. In fact, the cluster is already in that region 
at low temperature, as Fig. 7 shows. This is seen most 
clearly in the evolution of the diffusion behavior with 
temperature. In Fig. 8 we show the temperature varia- 
tion of the square root of the diffusion coefficient. While 
the two steps in the melting of the icosahedral isomer are 
clearly detected as slope changes at the corresponding 
temperatures, the value of \f~D for the disordered isomer 
increases with temperature in a smooth way, without any 
abrupt change. Thus, the opening of the available phase 
space proceeds in a gradual way and specific heat peaks 
are not detected. 

We have found a very different thermal behavior for 
two clusters that were classified in principle as disor- 
dered. Although just two examples are not enough to 
draw general conclusions, we believe that the thermal be- 
havior typical of amorphouslike sodium clusters is that 
found for Nag2, and that what is lacking is a clear defi- 
nition of what an amorphous cluster is. In line with the 
discussion of "atomic transfers" between shells advanced 
by Soler et al,cll we propose that a large orientational 
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disorder is not enough for a cluster to be classified as 
amorphouslike. Only when those atomic transfers are 
maximal, in such a way that no local regions with low 
atomic density exist, the cluster can be considered com- 
pletely amorphous. The existence of those regions, how- 
ever small they may be (as is the case of Nai42), promote 
the creation of appreciable free energy barriers in the po- 
tential energy surface, so sudden access to a substantial 
region of the available phase space is expected above a 
certain temperature, and peaks will appear in the tem- 
perature evolution of the specific heat. On the contrary, 
the absence of such low atomic density regions facilitates 
the diffusion of the atoms across the whole cluster volume 
right from the start of the heating process. As the liq- 
uidlike phase is established precisely when all the atoms 
in the cluster can diffuse across the cluster volume, we 
expect that no appreciable free energy barriers will be 
found in these cases, and no specific heat peaks will be 
detected. 

A few comments regarding the quality of the simu- 
lations and of the annealing runs are perhaps in order 
here. The orbital-free representation of the atomic in- 
teractions is much more efficient than the more accurate 
KS treatments, but is still substantially more expensive 
computationally than a simulation using phenomenolog- 
ical many-body potentials. Such potentials contain a 
number of parameters that are usually chosen by fitting 
some bulk and/or molecular properties. In contrast our 
model is free of external parameters, although there are 
approximations in the kinetic and exchange-correlation 
functionals. The orbital-free scheme accounts, albeit ap- 
proximately, for the effects of the detailed electronic dis- 
tribution on the total energy and the forces on the ions. 
We feel that this is particularly important in metallic 
clusters for which a large proportion of atoms are on the 
surface and experience a very different electronic envi- 
ronment than an atom in the interior. Furthermore, the 
adjustment of the electronic structure and consequently 
the energy and forces to rearrangements of the ions is also 
taken into account. But the price to be paid for the more 
accurate description of the interactions is a less complete 
statistical sampling of the phase space. The simulation 
times are substantially shorter than those that can be 
achieved in phenomenological simulations. The cooling 
rate employed in the annealing runs is also faster than 
those that can be achieved by using empirical potentials. 
Nevertheless, we expect that the locations of the several 
transitions are reliable, because all the indicators we have 
used, both thermal and structural ones, are in essential 
agreement on the temperature at which the transitions 
start. On the other hand, the disordered isomers found in 
different annealing runs did not show substantial struc- 
tural or energetic differences with respect to those studied 
here. 



IV. SUMMARY 

The melting-like transition in disordered Nai42 and 
Nag2 has been investigated by applying an orbital-free, 
density-functional molecular dynamics method. The 
computational effort which is required is modest in com- 
parison with the traditional Car-Parrinello Molecular 
Dynamics technique based on Kohn-Sham orbitals, that 
would be very costly for clusters of this size. This saving 
allows the study of large clusters. 

A disordered isomer of Nai42 melts in two steps as 
evidenced by the thermal indicators. The transition 
at T* s w 130 K is best described as intrashell melt- 
ing. This is followed at T* s w 270 K by homogeneous 
melting. Melting is found to depend on the starting 
low-temperature isomer. Specifically, for an icosahedral 
Nai42 isomer, the analysis of thermal, macroscopic prop- 
erties places those two stages much closer in temperature, 
240 K and 270 K respectively.Q Nevertheless, isomeriza- 
tion transitions are observed in the icosahedral isomer 
at a temperature as low as Tf IS «130 K. These isomer- 
izations involve the motion of the five atomic vacancies 
in the cluster surface, preserve the icosahedral structure 
and do not give rise to any pronounced feature in the 
caloric curve. In the disordered isomer there is not a 
separate isomerization regime (something rather evident 
because there is not an underlying ordered structure in 
each shell), but there is a melting-in-steps process, due to 
the distribution of the atoms in different shells. Thus, the 
melting of both isomers is not as different as suggested 
by the thermal indicators. An icosahedral isomer of Nag2 
melts also in a similar way: there is a minor peak in C v at 
rpico ~i30K indicates surface melting, and a main, 
homogeneous melting peak at T^° sa240 K. The lower 
value of T* co , as compared to Nai42, is due to the larger 
number of surface vacancies, and the melting-in-step pro- 
cess is due to the atomic shell structure. The melting of 
disordered Nag2 proceeds instead gradually, and spreads 
over a very wide temperature interval. The thermal indi- 
cators as the caloric curve or the specific heat do not show 
any indications of an abrupt transition, which suggests 
that the phase space available to the cluster does not in- 
creases suddenly at any given temperature. The square 
root of the diffusion coefficient increases with tempera- 
ture in a smooth way, in contrast to the step diffusive 
behavior of icosahedral Nag2 . It has been suggested that 
the absence of any abrupt transition is closely related to 
the absence of any shell structure in the radial atomic 
density distribution, which should be considered a neces- 
sary condition for a cluster to be considered completely 
amorphous. In this, sense, only the disordered isomer 
of Nag2 can be considered rigorously amorphous, while 
the disordered Nai42 isomer should be considered just 
partially amorphous. 

In summary, we have found a number of structural 
properties that have an important effect on the melting 
properties of sodium clusters. A melting in steps pro- 
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cess is to be expected in almost all clusters where a clear 
distribution of the atoms in radial shells exists, as is the 
case of both isomers of Nai42 and of the icosahedral iso- 
mer of Nag2- In those cases, intrashell diffusive motion 
starts at a temperature Ti ntra , lower than the temper- 
ature Tmter at which intershell diffusive motion begins 
to be important. The difference Tj n t er -Tj„t ro is small if 
the orientational order of the shells is large (for example, 
icosahedral order) and the number of vacancies in each 
shell is small: this occurs for icosahedral Nai42, with just 
five surface vacancies; A limit case is icosahedral Nass. 
With two complete atomic shells, intrashell motions are 
as difficult as ifttershell motions and the two transitions 
merge into one.0 When one shell contains a large number 
of vacancies, the two temperatures are well separated no 
matter how high the orientational order is: this is exem- 
plified by the case of icosahedral Nag2. Also, if the shells 
have a high orientational disorder, the two transitions 
are well separated in temperature no matter how close 
we are from an icosahedral shell closing: an example is 
the disordered Nai42 isomer. Finally, a gradual melting 
process without any abrupt transition is to be expected 
for all those clusters which have both orientational and 
radial disorder, that is for amorphouslike clusters: this is 
the case of amorphous Na 92 . 
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Captions of Figures. 

Figure 1 Structures of the low temperature (a) amor- 
phous Nai42, (b) icosahedral Nai42, (c) amorphous Nag2 
and (d) icosahedral Nag2 isomers. 

Figure 2 Short-time averaged distances < r,;(i) > sta 
between each atom and the center of mass in Nai42, as 
functions of time for (a) the icosahedral isomer at T= 30 
K and (b) the amorphous isomer at T= 47 K. 

Figure 3 Caloric and specific heat curves of Nai42, 
taking the internal cluster temperature as the indepen- 
dent variable. The deviation around the mean tempera- 
ture is smaller than the size of the circles. 

Figure 4 Short-time averaged distances < r,(t) > s t a 
between each atom and the center of mass in amorphous 
Nai42 as functions of time at T= 160 K. The bold lines 
follow the evolution of a particular atom in the surface 
shell and another in the outermost core shell. 

Figure 5 Time averaged radial atomic density distri- 
bution of the amorphous isomer of Nai42, at some repre- 
sentative temperatures. 

Figure 6 Caloric and specific heat curves of Na92 , tak- 
ing the internal cluster temperature as the independent 
variable. The deviation around the mean temperature is 
smaller than the size of the circles. 



Figure 7 Time averaged radial atomic density distri- 
bution of the amorphous isomer of Nag2 , at some repre- 
sentative temperatures. 

Figure 8 Square root of the diffusion coefficient as a 
function of temperature for the icosahedral and amor- 
phous isomers of Nag2 . 
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